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Abstract 

We have searched for the effective FCNC decays b — > sl + l~ using an inclusive 
method. We set upper limits on the branching ratios B(b — > se + e~) < 5.7 x 
10~ 5 , B(b -» sh+ij,-) < 5.8 x 10- 5 , and B(b -» se ± / u T ) < 2.2 x 10~ 5 (at 
90 % C.L.). Combining the di-electron and di-muon decay modes we find: 
B{b -» s/+r) < 4.2 x 10- 5 (at 90 % C.L.). 
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Flavor Changing Neutral Currents (FCNC) are forbidden to first order in the Standard 
Model. Second order loop diagrams, known as penguin and box diagrams, can generate 
effective FCNC which lead to b — > s transitions. These processes are of considerable interest 
because they are sensitive to Vt s , the Cabibbo-Kobayashi-Maskawa matrix element which 
will be very difficult to measure in direct decays of the top quark. These processes are also 
sensitive to non-Standard Model physics JT|, since charged Higgs bosons, new gauge bosons 
or supersymmetric particles can contribute via additional loop diagrams. 

The electromagnetic penguin decay b — > 57 was first observed by CLEO in the exclusive 
mode B -> K*-y with B(B -> K*j) = (4.2 ± 0.8 ± 0.6) x 10" 5 @. The inclusive rate 
for the decay B -> X s 7 was measured to be B(b -> S7) = (2.32 ± 0.57 ± 0.35) x f 0~ 4 
@. Inclusively measured rates are more interesting because they can be directly related to 
underlying quark transitions without large theoretical uncertainties in formation probabilities 
for specific hadronic final states. The measured inclusive b — > 37 rate is consistent with 
Standard Model calculations. 

The b — > sl + l~ decay rate is expected in the Standard Model to be nearly two orders of 
magnitude lower than the rate for b — > 37 decays. Nevertheless, the b — > sl + l~ process has 
received considerable attention since it offers a deeper insight into the effective Hamiltonian 
describing FCNC processes in B decays f|. While b — > 57 is only sensitive to the absolute 
value of the CV Wilson coefficient in the effective Hamiltonian, b — > sl + l~ is also sensitive to 
the sign of C7 and to the Cg and C10 coefficients, where the relative contributions vary with 
mass. These three coefficients are related to three different processes contributing to 
b — ► sl + l~: electromagnetic and electroweak penguins, and a box diagram. Processes beyond 
the Standard Model can alter both the magnitude and the sign of the Wilson coefficients. 
The higher-order QCD corrections for b — > sl + l~ are smaller than for the electromagnetic 
penguin and have been calculated in next-to- leading order 

Several experiments (UA1 0, CLEO 0, and CDF |J) have searched for the exclusive 
decays B — > Kl + l~ and B — > K*l + l~ and set upper limits at the level of (1 — 2) x 10~ 5 at 
90% confidence level (C.L.). Combining electron and muon modes, the previous generation 
of the CLEO experiment set an inclusive limit: B(b -> sl + l~) < 1.2 x 10" 3 (90 % C.L.) 0. 
The UAl experiment J7| searched for inclusive b — > sfi + fi~ at the endpoint of the dilepton 
mass distribution (M(fi + yr) > 3.9 GeV) which comprises about a tenth of the total rate. 
Extrapolating to the full phase space UAl claims a limit of < 5 x 10 -5 (90 % C.L.). However, 
a simulation of the UAl acceptance shows that UAl overestimated their efficiency by at least 



a factor of three [HT 



In this Letter, we present results of the search for inclusive b s/i + yU~ , b — > se + e~ and 
b — > se ± fi T . The latter decay violates conservation of electron and muon lepton numbers 
and thus can only originate from processes beyond the Standard Model. The data were 
obtained with the CLEO II detector at the Cornell Electron Storage Ring. A sample with 
an integrated luminosity of 3.1 fb _1 was collected on the T(4S') resonance. This sample 
contains (3.30 ± 0.06) x 10 6 produced BB pairs. For background subtraction we also use 
1.6 fb _1 of data collected just below the T(4S*). CLEO II is a general purpose solenoidal 
spectrometer described in detail in Ref. fl2fl . 

The data selection method is very similar to the reconstruction method presented in our 
previous measurement of the b — > S7 rate 0, with the 7 candidate replaced by a lepton 
pair. We select events that pass general hadronic event criteria based on charged track 
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multiplicity, visible energy and location of the event vertex. The highest energy pair of 
oppositely charged leptons is then selected. Electron candidates are required to have an 
energy deposition in the calorimeter nearly equal to the measured momentum, and to have 
a specific ionization (dE/dx) in the drift chamber consistent with that expected for an 
electron. Muon candidates are identified as charged tracks with matching muon-detector 
hits at absorber depths of at least three nuclear interaction lengths. In the channel, 
one muon is required to penetrate at least five interaction lengths. We then look for a 
combination of hadronic particles, denoted X s , with a kaon candidate and 0-4 pions, which 
together with the selected lepton pair satisfy energy-momentum constraints for the B decay 
hypothesis, B — > X s l + l~ . To quantify consistency with this hypothesis we use: 

2 [Mb — 5.279 \ 2 / Eb — Eb eam 



Xb = 

V cr M J \ a E 

where Mb = \J El eam — P§, Eb-, Pb are the measured energy and momentum of the B 
candidate, and <tm, &e are experimental errors on M B and E B estimated from the detector 
resolution and beam energy spread. The kaon candidate is a charged track with dE/dx 
and time of flight (ToF) consistent with the kaon hypothesis, or a K s — > -k + -k^ candidate 
identified by a displaced vertex and invariant mass cut. A pion candidate is a charged track 
with dE/dx and ToF consistent with the pion hypothesis, or a 7r° — > 77 candidate. At 
most one 7r° is allowed in the X s combination. In each event we pick the combination that 
minimizes overall \ 2 1 which includes Xb together with contributions from dE/dx, ToF, and 
K s and 7r° mass deviations, where relevant. 

To suppress continuum background we require the event to have H 2 /H Q < 0.45, where Hi 



are Fox- Wolfram moments ||13|| . We also require | cos 9 tt \ < 0.8, where 9 U is the angle between 
the thrust axis of the candidate B and the thrust axis of the rest of the event. To suppress 
BB background we require the mass of the X s system to be less than 1.8 GeV. The dominant 
BB background comes from two semileptonic decays of B or D mesons, which produce the 
lepton pair with two undetected neutrinos. Since most signal events are expected to have 
zero or one neutrino, we also require the mass of the undetected system in the event to be 
less than 1.5 GeV. By excluding the mode with X s = Ktt+tt^tt , we reduce the expected 
BB background by an additional 21% while reducing the signal efficiency by only 6%. 

Fig. |I] shows the dilepton mass, M(l + l~), for the events which pass the cuts previously 
described and the B consistency requirement, Xb < 6, in the on- and off-resonance data 
samples. Unlike the b — > 57 analysis, the continuum background is small. The peaks at the 
i/j and ip' masses that are observed in the on-resonance data are due to the long distance 
interactions, b —> cW~ , W~ — > cs, with the cc pair hadronizing into a if)^ which then decays 
to a lepton pair. Using cuts on M(X S ) to identify K and K*, we measure the branching 



ratios for B — > ipC'K^*' and obtain results consistent with a recent CLEO publication [|T4" | 
For further analysis, we exclude events with M(l + l~) near the ip and ip' masses (±0.1 GeV 
for —0.3, +0.1 GeV for e + e~, no cut for e ± fi T ). The exclusion region is wider in the 

e + e _ channel because of the radiative tail. After these cuts and continuum subtraction, we 
observe 10 ± 5 X s e + e~ , 12 ± 6 X s/ u + /i~, and 18 ± 8 V s e ± /i =F events in the data. From the 
Monte Carlo simulation of generic BB events we expect 9 ± 1, 16 ± 2 and 39 ± 3 (statistical 
errors only) background events, respectively. Therefore, no evidence for signal is found in 
the data and we proceed to set limits on these decay rates. 
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FIG. 1. M{l + l ) distributions for the on- (top) and off-resonance (bottom) data with the 
x\ < 6 cut. The scaling factor between the off- and on-resonance data is 1.9. 

To avoid systematics related to absolute normalization of the BB Monte Carlo, instead 
of counting events after the Xb < 6 cut, we loosen this cut to 30 and fit the observed Xb 
distributions in the on- and off-resonance data. We allow for signal contribution, as well as 
BB and continuum backgrounds. The signal is expected to peak sharply at zero, whereas the 
backgrounds have natter distributions (see Fig. (a)). The shapes of all these contributions 
are fixed from the Monte Carlo simulation, while the normalizations are allowed to float. The 
relative normalization for continuum background between the on- and off-resonance data is 
fixed to the known ratio of integral luminosities and cross-sections. The fitted number of 
X s e + e~ , X s fi + fi~ , and X s e ± /i =F events is 7 ± 7, 1 ± 7, and —18 ± 10, respectively. As an 
example, the fit to the X s e + e~ data is displayed in Fig. ^| (b)-(c). 

To calculate the signal efficiency and to predict the x% signal distribution we gener- 
ated b —>■ sl + l~ Monte Carlo events. The parton level distributions for b — ► se + e~ and 
b — ► sjj^yT are predicted from the effective Hamiltonian containing Standard Model con- 
tributions. The next-to-leading-order calculations were used |H|. At present, the effect of 
gluon bremsstrahlung on the outgoing s— quark is only partially included in the theoretical 
calculations. After our ip and if)' veto cuts, the long distance interactions are expected to 
constructively interfere with the short distance contributions. Estimates of these interference 
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effects are model dependent. The most recent calculation predicts modifications of the short 
distance rate by only about 2% |ff5] compared to 20% predicted by some earlier simplified 



models | I6| . We neglect long distance interactions in our Monte Carlo. Since no theoreti- 
cal calculations for the non-Standard Model decay b — > se ± /i =F exist, we use a phase space 
model for these decays. To account for Fermi motion of the b quark inside the B meson 
we have used the spectator model by Ali et al. fl7f . The particle content of the X s system 



was modeled with the conventional method of quark hadronization from JETSET |Tj|. For 
better accuracy of the simulations, when M(X S ) is in the K or K* mass region, the event is 
regenerated according to the theoretical predictions for the exclusive B — » K^*H + l~ decays 
by Greub et al. |0|. The estimated efficiencies are 5.2%, 4.5% and 7.3% for e + e~, and 
e ± /i =F modes respectively. 

To estimate the systematic error due to the uncertainty in the Xb signal and background 
shapes, we divide the Monte Carlo sample into low and high multiplicity channels in the 
manner which produces the largest shape variation. This shape variation changes the upper 
limits by 9%, 19% and 20% for the e + e~, and e ± n^ channels respectively. Uncertainties 

in detector simulation are dominated by charged tracking systematics and are estimated to 
be 14%. Variations of the spectator model parameters || result in changes of the selection 
efficiency by (12±4)%, (30±4)%, and (11 ±4)%, respectively. The larger uncertainty in the 
channel is the result of the lack of muon identification for P M < 1 GeV/c. Uncertainty 
in the modeling of the hadronization of the X s system gives a contribution of 9%. Adding 
these sources of systematic errors in quadrature, we estimate the total systematic errors to 
be 22%, 39% and 28%, respectively. 

Using a Gaussian likelihood integrated over positive signal values, we find upper limits 
using statistical errors only. We then loosen these limits by one unit of systematic uncer- 
tainty. 

The final results are B{b — > se + e~) < 5.7 x 10~ 5 , B(b — > s/i+//~) < 5.8 x 10~ 5 , and 
B{b — > se ± /x T ) < 2.2 x 10~ 5 . The results are consistent with the Standard Model predictions 
03, (0.8 ± 0.2) x 10^ 5 , (0.6 ± 0.1) x 10^ 5 , and 0, respectively. Combining the e+e" and 
results, we also set a limit on the rate averaged over lepton flavors, B(b — > sl + l~) = 
(B{b -> se+e-) + B(b -> s/i+/i~))/2 < 4.2 x 10" 5 (90 % C.L.). 

The limit on B(b — > se + e~) is more than an order of magnitude more restrictive than 
the previous limits. The limit on B(b — > sfi + [i~) is also significantly tighter than the UAl 
limit after correcting for the efficiency problem (see discussion above). Furthermore, in 
contrast with the UAl measurement, the present analysis is sensitive to a much wider range 
of M(/ + / _ ). Therefore, our extrapolation to the full phase space is more reliable and we are 
sensitive to a broader range of processes beyond the Standard Model. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the National Science Foun- 
dation, the U.S. Department of Energy, the Heisenberg Foundation, the Alexander von 
Humboldt Stiftung, the Natural Sciences and Engineering Research Council of Canada, and 
the A. P. Sloan Foundation. 
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FIG. 2. Xb distributions for X s e + e~ data, (a) The difference between the expected distribution 
for the signal (solid histogram) and BB background (dashed histogram). Both distributions are 
normalized to the same area, (b) The fit to the on-resonance data (points with error bars). The 
sum of all fitted contributions is indicated by a solid line. The fitted background contribution (BB 
plus continuum) is indicated by a dashed line. The estimated continuum background, indicated by 
a dotted line, is simultaneously constrained to the off-resonance data (c). 
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